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We present a method to determine the nodal structure of the energy gap of unconventional su- 
perconductors such as high T c materials. We show how nonlinear electrodynamics phenomena in the 
Meissner regime, arising from the presence of lines on the Fermi surface where the superconducting 
energy gap is very small or zero, can be used to perform "node spectroscopy" , that is, as a sensitive 
bulk probe to locate the angular position of those lines. In calculating the nonlinear supercurrent 
response, we include the effects of orthorhombic distortion and a — b plane anisotropy. Analytic 
results presented demonstrate a systematic way to experimentally distinguish order parameters of 
different symmetries, including cases with mixed symmetry (for example, d + s and s + id). We con- 
sider, as suggested by various experiments, order parameters with predominantly d-wave character, 
and describe how to determine the possible presence of other symmetries. The nonlinear magnetic 
moment displays a distinct behavior if nodes in the gap are absent but regions with small, finite, 
values of the energy gap exist. 
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I. INTRODUCTION 



In recent years a very significant effort, both experimentally and theoretically, has been made to determine the 
symmetry of the pairing stateou in high temperature superconductors (HTSC's). It is usually held that such an 
identification would provide significant clues towards a better understanding of the microscopic mechanism respon- 
sible for superconductivity. The symmetry properties of the pairing state itself, being connected to those of the 
superconducting energy gap, have important physical and technological implications. 

There is a large body of experimental results, supported by theoretical work, that is generally interacted as 
indicative of a pairing state at least predominantly of the d-wave typeutffl, (although there are experimentsEra which 
are difficult to interpret in that framework.) By this, it is meant a pairing state which has an order parameter 
(OP) which vanishes, or nearly so, along four lines on the Fermi surface (FS). The angle between such lines of nodes 
(or "quasinodes" ) is presumed to be near 7r/2, but, because of the, orthorhombic symmetry of HTSC materials (as 
exemplified by the strong anisotropy of superconducting propertiesEij in the a — b plane) this is unlikely to be precisely 
correct. Pairing states of different tetragonal symmetry, such as d- or s-wave, could nip- -Several proposals, mostly 
incorporating a predominantly c£-wave OP, have been set foeth to describe this mixingE2rt-X Some of them include 
OP's of different symmetry near the surface than in the bulkllj. Many of the best probes of the OP give information 
only on the pairing state within a few correlation lengths from the surface, which might differ from the bulk state. 
Thus, large uncertainties exist. I— . .— . 

The aim of this paper is to show that nonlinear electromagnetic phenomenaliH~E3 can provide a systematic way 
to distinguish OP's of different symmetries, or of mixed symmetry. We will show how the nonlinear electrodynamic 
effects in the Meissner regime, which are due to the presence of nodes or small minima in the superconducting gap, can 
be used to perform "node spectroscopy" , that is, to determine the angular position of any lines on the FS where the 
absolute value of the OP is very small or zero, and the value of the OP at these lines. Further, the bulk character of 
these effects (they extend over a length on the order of the penetration depth rather than the much smaller coherence 
length) means that these determinations would refer to the bulk, and can become very useful in avoiding material 
imperfections near the surface and other difficulties of the surface probes. 

We concentrate in this work on the nonlinear magnetic moment as a sensitive probe for the superconducting gap 
symmetry. In particular, we examine its transverse component (perpendicular to the applied field, which lies in 
the a — b plane) which is determined by the quasiparticle excitations. Other quantities, such as the penetration 
depth, might also be used for the same purpose, and will be referred to briefly The transverse moment, however, 
can be readily measured either directly or through the torque it produces, and, as shown in Ref. [l9[ its size for 
typical samples is large enough to be detected. At temperatures much below the critical temperature, the excitations 
responsible for the nonlinear Meissner effect are energetically favorable only near the nodes or quasinodes of the FS. 
The spectroscopy method described in this paper, provides an accurate positioning of the nodes in the energy gap 
and not just information about their presence. Similarly we show that even in the absence of nodes, the presence of 
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"quasinodes" would lead to several distinctive features in the angular and field dependence of the nonlinear magnetic 
moment. 

We will examine various pairing states that are compatible with the crystal symmetry, including leading candidates 
for the paiting state of HTSC materials, and excluding only the pure d-wave state which we have extensively discussed 
elsewheretj. In our notation for the pairing states, we shall follow the discussion from Refs. [l]j|. We start with the 
spin singlet pairing states in a square CuOi lattice as given in Table |[ For the YBCO family of cuprates, the 
orthorhorpbic crystal distortion amounts to a simple change in the length of the a and b crystallographic constants. 
Thereforeo the admissible mixing of OP's include combinations of d x z_ y 2 and s, or of s~ and d xy symmetries only. 
The nodes of the gap function are not necessarily at an angle 7r/4 from the cpzatallographic axes, and do not in 
general form an angle of tt/2. For the BSCCO family the situation is differentEH. The distortion occurs in such a 
way that the new orthorhombic axes form an angle of 7r/4 with the original tetragonal axes. The relevant mixing is 
then between d xy and s or between d x 2_ y 2 and s~ symmetries. The nodes are at the orthorhombic axes and, as in 
the pure <i-wave case, the directions along which they occur form angles of tt/2. In the following sections, we shall 
use these considerations to choose the gap functions to be studied. 

In Section II we discuss methods to solve the nonlinear Maxwell-Equations in the Meissner regime to compute the 
physical quantities of interest. The general procedure we use exploits the fact, proved in Ref. ^| (see particularly Fig. 
5 there) that the angular dependence of the transverse magnetic moment for a typical finite sample (a flat crystal) can 
be very accurately computed by assuming that the sample is infinite in the a — b plane, while the overall amplitude is 
underestimated somewhat by this assumption. Since, as we shall see, the nodal structure is directly reflected in the 
angular dependences, we can proceed for our purposes here by considering a simple "slab" geometry. This allows us 
to do the entire work analytically, treating the small nonlinear effects perturbatively. We illustrate our perturbation 
method in Section II by considering first a superconducting gap with mixed d + s symmetry. Various additional gap 
functions, with nodes, are studied in Section III. We show how to include the effects of orthorhombic distortion as they 
affect both normal (FS shape) and superconducting (penetration depth tensor) properties, in the calculation of the 
nonlinear effects. In Section IV we consider gaps without nodes in contrast with the corresponding results obtained 
in the previous sections. We propose several ways that would allow to experimentally distinguish nodes from small 
minima in the superconducting gap. In the last Section we give our conclusions and possible guidelines for the future 
work. 



II. METHODS 



A. Maxwell-London Electrodynamics 

We briefly review here the nonlinear Maxwell-London equations&S'il that need to be solved for superconduc- 
tors with various pairing states. For a d wave infinite ^superconductor, these nonlinear equations have been solved, 
perturbativelju at zero temperature, or numericallyOu when finite systems or temperatures have been considered. 

Here we are concerned with node spectroscopy, that is, with the effect that a pattern of nodes (or deep dips) in the 
angular dependence of the OP has on the dependence of the nonlinear Meissner effect OB-jbhe angle between applied 
field and sample orientation. For this purpose we can work at T = since it was shownO that at the temperatures 
T «T C and field ranges experiments are performed, the T = result iSpSufficiently accurate. Similarly, we can work 
in a simple "infinite slab" geometry (as defined below) since it is knowrfiia that finite size effects, while changing the 
amplitude of the nonlinear phenomena to some extent, do not affect their angular dependence. We will, therefore, 
develop a zero temperature perturbation methods, for a simple geometry, that can be used with a wide variety of 
pairing states. These methods ate,|6unontrivial extension of those used in Ref. [l8| for a pure d-wave state. 

We consider a superconductortlnia infinite in the a — b plane and having a thickness d in the c direction. We will 
assume unless otherwise stated that d is larger than any relevant penetration depth, (a "slab", not a film). We use 
a coordinate system with z axis parallel to the c crystallographic direction. We follow the notation introduced in 
Ref. [l9| unless otherwise indicated. The superfluid velocity v is defined as: 

v=^ + -A, (2.1) 
2 c 

where x is the phase of the superconducting OP, A the vector potential, and e the proton charge, (with Ti = ks = 1). 
The relation between v and H is given by the second London equation: 

Vxv=-H. (2.2) 

c 
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In the steady state the appropriate Maxwell equation is Ampere's law, V x H 
we obtain: 

_ _ 47re., . 

V X V X V = — 5-jlv). 



The relation between j and v is generally nonlinear and given by 



^j. Combining it with Eq. 

(2.3) 



j(v) 



eN f f d 2 S n( S )v / [(v / -v) + 2 / dl; /(£(£) + Vj • v)], 
Jfs Jo 



(2.4) 



where Nf is the total density of states at the Fermi level, n(s) the density of states at point s at the Fermi surface 
(FS), normalized to unity (J FS d 2 s n(s) — 1), v/(s) the s-dependent Fermi velocity, / the Fermi function, with 

E($) = (£ 2 + lA(s)! 2 ) 1 / 2 , T the absolute temperature and A(s) the OP. 

For the geometry considered and for the magnetic field H a , applied in the a — b plane, the fields have only x and y 
components, which depend only on the coordinate z. Eq. (2.3) then reduces to 



o zz v + — 5-j(v) 



0. 



where it is understood that Eq. (2.4) is substituted in the second term. The boundary conditions are: 

H = H a | z=± d. 



(2.5) 



(2.6) 



The first term in ( 2.4) is the usual relation j = — epv, w here p is the superfluid density tensor. At T = the nonlinear 
corrections (in v), described by the second term of Eq. (2.4), and due to the quasiparticle backflow can be written as: 



j„„(v) -2,X f j ./-',- /,(.s) V/ B. - V/ .v -|A ( .i|i[,v,.v:r - |A(*)| 2 ] 1/2 

n= 1,2,.. 



FS 



2eE^/« / #^ /n [(v /n • v) 2 - |A(^)| 2 ] 1/2 , 



(2.7) 



where we have replaced integration over the whole FS by integration over only the nodal (or quasinodal) regions fi n , 
each region enumerated by the index n, where the gap vanishes (or dips to a value much smaller than its maximum) 
so that a backflow current due to so calledEJ "jets" of quasiparticle excitations is allowed. Quantities at the node n 
are denoted by an index "n". For an anisotropic FS surface some care has to be taken: at the nodal position s n , 
the density of states, n(s„), is not necessarily equal to the average value, rii S0 . At each node we therefore define 
the appropriate weighted density of states Nf n = Nfn(s n )/rii S0 . Regions contributing to the backflow current are 
determined by 



v/ • v < 0, 



(2.8) 



where <f> n is angle with respect to the closest node, n, in the energy gap, and can be approximated by v/ n , its value 
at that node. As an example, for a pure d-wave (d x 2_ y 2) OP, A(<fi) = A^ sin(20), where <fi is measured with respect 
to a node, Eq. (2.7) can be integrated with f2„ given from Eq. (2.8) by \cf> n \ < 



( V fr, 



|A d sin20 c 



4A 



2^,2 
cn ' 



cn , where cj) cn is determined from: 

(2.9) 



and <f> cn <C 1, which implies that the nodal value is an accurate approximation on VL n , It follows from these consid- 
erations that the nonlinear effects depend on the local values of FS quantities at the nodes. The overall shape of the 
FS enters (as we shall see in more detail below) very indirectly, through the eigenvalues of the London penetration 
depth tensor. 



B. Nonlinear Currents 



To perform the calculation of the angular dependence of the nonlinear electrodynami c q uantities we must first 
develop a general method to compute the current, including the nonlinear terms from Eq. (2.7), for the different OP's 
under consideration. The resulting expressions lead, after substitution into (2.5), to nonlinear differential equations 
that can be solved perturbatively. It is pedagogically convenient to illustrate both steps by carrying them out in 
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terms of a mixed d + s-wave gap function and an isotropic FS. It is then rather straightforward to extend the results 
to the other forms of the OP and to include the effects of the a — b plane anisotropy. 
The specific form of a d + s OP that we consider in this paper is: 



(<j>) = A s + A rf sin20, (2.10) 



where all quantities are real, <fi is measured from the X axis (Fig. [I]), and we take |A S | <C |A<j|, which ensures that 
there are nodes in the gap,£imction. The experimental constraints on the smallness of A, / A^ have been discussed, for 
example, by Annett et al.Ero. At A s = we will recover the known perturbative resultsES for that limit. The relevance 
of this form of the gap to YBCO, wh en th e effects of orthorhombicity are taken into account, has been discussed 



above. For A s 7^ 0, the nodes of Eq. ( 2.10 ) are no longer separated by an angle of tt/2 since they are shifted by an 
angle ±v (see Fig. |l|), 

from the orthogonal axes, which we denote by X and Y, and which will remain defined throughout the paper as being 
along the nodal directions of a pure d-wave gap, that is, at angles of ir/4 with the crystallographic axes of YBCO. 
These axes, together with other details, are shown in Fig. [lj We introduce also nonorthogonal axes along the nodal 
directions of the d + s gap function with unit vectors x n (n=l,2,3,4 as indicated in the Figure): 

x\ = -x 3 = cX - sY, x 2 = -x 4 = -sX + cY, (2.12) 



where c = cos is and s = siniA The angles <p n , defined below Eq. ( |2.8| ), measure angular shifts from the nodes of 
Ad+s, for example, <pi = 4> + v. Close to the nodal region, the gap has the form: 

|A rf+s (0„)| m |/iA e// n | , n = l,2,... (2.13) 

where the effective amplitude of the gap function A e // = (A^ — A;;) 1 / 2 is slightly reduced from the pure d-wecve limit, 



and n = 2 in this case. The allowed region for quasiparticle excitations from Eq. (2.8) is §\ < <j)^, n where 

We will investigate the angular dependence of physical quantities, when the field is applied along different directions 
with respect to nodes in the gap function. We defin e, as in R ef. |l9| , ip to be the angle of H a with the +Y direction. 
This angle is also shown in Fig. 0. From Eqs. (§75]) and (fl| we see that as ip changes, various nodes contribute to 



j<7P- Without a — b plane FS anisotropy the integration yields (as a special case from Appendix ^) 

1 v 3 f 

jgpi = T; eN f a \ e u{f)tr x + e 2 i{^)v x VY + e 3i (ip)v^], i = X,Y, (2.15) 

2 (lAeff 

where the coefficients ei j 2,3i(V') are defined for various ranges of ip: 

eix = e 3Y = c 3 - s 3 , e 2 x = 2e 3X = 2e 1Y = e 2 y = ~2cs(c - s), ip € [v, ~ - v] (2.16a) 

eix = -e 3 y = -(c 3 + s 3 ), e 2X = -2e 3X = 2e 1Y = -e 2Y = 2cs(c + s), ip G - v, ir + v). (2.16b) 

We can simply verify that for a pure d wave case (c = 1, s = 0) these expressions correctly reduceEl to: 

v x ,y\v x ,y\ 

J qP x,Y = ep ab — , (2.17) 

where p a b = is th e superfluid density, related to the penetration depth by \~ b 2 = ^§-p a b, and the characteristic 

velocity in that limit is v c = Ai. I n general, we define 

v c ee A eff /v f , (2.18) 
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and a dimensionless velocity 



m = , i — X,Y. 

jJLVc 



(2.19) 



We also define the dimensionless field: 



Ha 
Hn 



(2.20) 



where we introduce the characteristic magnetic field as: 



H 



cA 



eff 



eXabVf 



In the case of pure <i-wave gap, Eq. (2.21) reducesEl to Hq = ^2\° b ^ - 
coherence length is defined as £ a b = Then, recalling (2.15), we can write (2.5) in dimensionless form as: 



(2.21) 



, where 0o is the flux quantum and the in-plane 



dzzUi ~Ui + [eii(ip)u 2 x + e 2 i(ip)uxu Y + e 3i (V>)uy] = 0, 



X,Y. 



(2.22) 



Here we have introduced the dimensionless coordinate Z = z/X ao , where X a b is the penetratipa depth at any of the four 
equivalent nodes. In the pure <i-wave limit the above equations again reduce to the knownEEl result. In the geometry 
considered in this paper j qp i and Ui have definite parity with respect to the Z coordinate and it is sufficient to solve 
the boundary value problem for half of the slab (Z > 0). The boundary conditions at the surface Z = Z. 
from Eq. (2.2) are: 



2A Q 



dzu x \z=z s = h 



cost/' 



dzu Y \z=z, = h 



sin-0 



(2.23) 



and odd parity of u requires: 



ux,y 



0| 



(2.24) 



C. Perturbation method 



In the Meissner regime where h, u <C 1, Eq. (2.22), with the boundary conditions ( 2_23| ) and ( ^.24 ), can be 
solved perturbatively. We write Ui(Z) as a sum of linear and quadratic parts in the smalEJ parameter h, Ui(Z) = 
Li(Z) + Ni(Z). The linear part then satisfies: 



d zz L i (Z)-L i (Z)=0, 



with Li satisfying the boundary conditions ( 2.23 ), ( 2.24 ). The solution is: 

L l (Z) = ai sinh(Z) , i = X, Y 



and from Eq. (2.23) 



a x = h- 



cosip 



[i cosh(Z s ) ' 

The nonlinear, perturbation part, in h, Ni, satisfies : 



ay = h 



sin-0 



/xcosh(Z s 



dzzNi(Z) - Ni{Z) + aj sinh 2 (Z) = 0, 



af = [eu(ip)a x + e 2l {^)a x a Y + e 3i (ip)al 



X,Y 



at Z = Z s from Eq. (2.23) the boundary conditions are: 

d z N l \ z =z s = 0. 



(2.25) 
(2.26) 

(2.27) 

(2.28a) 
(2.28b) 

(2.29) 
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The local magnetic field at Z = Z s is equal to the applied field. At Z = 0, we have from Eq. (2.24) that A^(0) = 0. 
The complete solution to Eq. j2.2Sf) is a sum of the homogeneous solution, Nh% and a particular solution, N p i\ 



Ni(Z) = An sinh(Z) + A 2l cosh(Z) + P, cosh(2Z) + R t = N hl {Z) + N vi (Z), 



(2.30a) 



A u = \a\ (sinh(Z s ) - tanh(Z,)), A 2t = -\a 2 , P i = \n i = - a ±, i = X,Y. 
3 3 3 o 



(2.30b) 



where the coefficients in iV^i, An, A 2i , are obtained using Eq. ( |2~24] ) and ( gTgg ). The magnetic field can be calculated 
from the field u through Eq. (2.2): 



H X (Z) = - H Q [hf L {Z) sin ijj + — (e 1Y cos 2 ip + e 2Y cos-0 sin0 + e 3Y sin 2 ip)f N (Z)}, 



(2.31a) 



H Y (Z) = +H Q [hf L {Z) cos?P + —(e ix cos 2 i/j + e 2X cosiPsm^ + e 3X sin 2 -iJj)f N {Z)} (2.31b) 

where 

h(Z) = ^P-, (2.31c) 
cosh(Z s ) 

f N (Z) = ^ [ 2(sinh(Z s ) - tanh(Z s )) cosh(Z) + 2sinh(Z) - sinh(2Z)]. 

3 cosh (Z s ) 

The above results are valid for any value of d (not just (d 3> A)), and we can at once write expressions for measurable 
quantities. The components of the penetration depth along the two orthogonal directions defined by the axes can be 
obtained as: 

1 = ^-\d z H x \z=z., . * = ^-\d z H Y \z=z.. (2.32) 



Aj>f (H a ) H a ' Ay(H a ) H a < 

For a slab, Z s 3> 1 (rf > A), a pure <i-wave gap function (v — 0) and field applied along the node or antinode, one can 
use these and similar expressions to recover the results for l/A(H a ) from Refs. |l7|]lS| . 

We next consider the magnetic moment, in particular its component transv erse t o the applied magnetic field, which 
in this case is purely nonlinear for a symmetric samplec.il . In previous workMlE3 we have shown that the magnetic 
moment (for any sample geometry) can be calculated using only the surface values of the fields. For the geometry 
considered here, the magnetic moment, defined in general as: 

m=iy"drrxj(v), (2.33) 

can be expressed as: 

SdH axy Sc . . . /r>o^\ 

m *,y = 4^ T 2TO %a ^/^' ( ) 

where x and y are orthogonal axes fixed in space, S is the.slab area, and we have used that v i s odd in zil. For a 
slab, (d ^> A ao ), which is the case of experimental interestO, the expressions for N X ^ Y from Eq. ( 2.30 ) simplify: 

h 2 

Ni(Z s ) = — —(eu cos 2 if) + e 2l cos -0 sin ^ + e 3i : sin 2 ijj) i = X,Y (2.35) 
and we obtain the transverse magnetic moment, mi = m x cosip + m Y simp: 



SX ab H 2 
3[i 2 ir Ho 

+ cos ip sin ip((e ix - e 2Y ) cosip + (e 2X - e 3Y ) sin if))]. 



"U WO = sin3 ^ - eiv cos 3 (2.36) 



The most obvious consequence of the s-wave admixture is that the periodicity of m±(t/)) is no longer strictly 7r/2, but 
it, although the 7r/2 Fourier component remains larger for small v. Writing explicitly ei,2,3i we have: 



G 



Tn±(ip) = -77T~ ~rr( c ~~ s)(cosip — sin ip)[cs + (1 + Acs) cos ip sin ip] , )/) e [i/, ^ — i/] (2.37a) 

127T Hq 2 

mx(V') = ~ ^ n ab tt( c + s)(cos^/> + sin ip)[cs + (1 — 4cs) cos f/ 1 sin ^] , ip € [— — z/, n + v\. (2.37b) 
127T Hq 2 

We can readily extend this solution for m±(ip) to the remaining ip tange by symmetry. In the limit v — ► 0, corre- 
sponding to a pure <i-wave, we recover the known result for this limito. 

The above expressions are in terms of the angle ip with a line of nodes. In an experimental situation the position of 
the line of nodes is not a priori known, since it depends on the amount of A s admixture, which is to be determined. 
All one has to do is to plot mi as a function of angle, with arbitrary origin, and compare the shape of the curve, up 
to an horizontal translation, with the plots presented below. 

We now show in Fig. || the angular dependence of m±(ip) from Eq. (2.37) for ip S [— f , %-]■ This interval is 



a complete period for the plotted quantity. The solid line represents the pure d-wave result, with its maximum 
normalized to unity. The actual unnormalizcd value of its amplitude as a function of applied field and sample 



parameters can be read off Eq. ( 2.37 ). It is discussed in detail, with inclusion of finite size effects, in Refs. |19|j20| , 
where typical values of Hq, h, S, and other sample parameters, arc given. The other lines depict m±(ip) for a d+s gap. 
They are normalized by the same factor as the solid line, so that the relative effects of introducing an s component 



can be gauged. These curves, which are (see Eq. ( [2.37 )) independent of the applied field, are labeled by the value of 



the a ngle v, the shift of the nodal directions from the ci- wave nodes. The relation between v and A s /Ad is given by 



Eq. ( 2.11 ). If we take v < (d — s gap) the lines are simply shifted by n/2. In Fig. |2j we see that when the nodal 
directions are not orthogonal, and hence the period of m± doubles from its d-wave value, the formerly equivalent 
values of the maxima and minima of |rnx(^)| take increasingly different values. One of the maxima becomes now the 
primary maximum, and its value increases with while the value of the other maximum decreases with \v\. The 
overall increase in the maximum nonlinear effect occurs largely because the contributions from the two quasiparticle 
"jets" at each node tend to cancel each other to a lesser extent when the angle between them is larger than tt/2, that 
is, near the primary maximum. From the symmetry of the d + s gap function, for any v (and isotropic FS), when the 
field is applied along antinodes, v || j and m±(ip) = when ip = 2{n + 1)tt/4, n = 0, 1, 2, .. (for <i-wave this is also true 
for the field applied along the nodes). 

For experimental analysis, it may be convenient to examine the main Fourier components of the periodic m±(ip) 
curve, in particular its part quadratic in the applied field H a . As one can see from Fig. ||, the transverse moment has 
odd parity with respect to ip if the angle is measured from a line of antinodes. With that choice, only the coefficients 
of the sine terms in the Fourie r analysis are nonvanishing. The ones corresponding to 7r and 7r/2 periodicities are 



easily obtained from Eq. ( 2.37 ): 

SX H^ 111 

ni2 = — -[cs(l H — cos 2v) H sin 2v sin 6v], (2.38a) 

37r Hq 5 12 20 

SXabH^l 2c 2 s 2_ 2 cog4l/+ 1 cos6u+ * ssin62 ,] i (2.38b) 



7714 



3tt ff 4 15 28 7 



where mi and and 777,4 are the 7r and ^ Fourier components of m±(ip), respectively. When Fourier components are 
computed with respect to an arbitrary angular origin, one will in general measure both sine and cosine co effici ents, 
but the square root of the sum of their squares should equal the absolute value of the results given in Eq. ( 2.38| ). 



In Fig. |^ we plot the Fourier components 7712, 7774, as a function of v and at a constant field, normalized so that 7774 
is unity for d-wave. The solid line represents 7774. We see that despite the overall period having doubled, this Fourier 
component is only slightly reduced for the range of v we considera. The n component, 7712, is depicted by the dotted 
line and it increases monotonically with v. 

In the following sections we will focus on the nonlinear transverse magnetic moment for different forms of gap func- 
tions. Other quantities could be obtained as shown above. We will see that the generalization of the methods discussed 
in this Section for the evaluation of the currents and for the perturbation values of the fields is straightforward. 

III. GAP FUNCTIONS WITH NODES 

In this section we consider additional relevant order parameters with nodes. We will include the details of the 
derivations in so far as they differ from the example discussed above, but the emphasis will be in the results for the 
angular dependence of the transverse moment. 
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A. Anisotropic d-wave states 



The first examples we consider are those of a pure ci-wave order parameter, with nodes at n/2 angles, but where 
the FS in the a — b plane is anisotropic, so that the penetration depths are not equivalent. This is a possible state 
for both the YBCO and BSCCO families. For the latter case the situation is easier since the Fermi velocities at the 
nodes are still aligned with the X or Y orthogonal axes. We therefore concentrate on the more complicated first case, 
where the Fermi velocities at the nodes are not aligned with these axes, but are equal in magnitude. 

1. YBCO type orthorhombicity 

For YBCO, although the oi±horhombic distortion of the crystal lattice is small, there is a quite significant in-plane 
penetration depth anisotropyJiJ (A 2 /A 2 ^2 — 6). Even for a pure d-wave gap function, the vectors V/„ no longer lie 
along the nodal directions ±X, ±Y, but are shifted by an angle ±1/4- We define va as the angle between the +Y 
axis and vp. In Fig. [| we show relevant quantities for an anisotropic FS and the more general case of a d + s gap. 
The pure d-wave limit, considered in this subsection, corresponds to setting v = in the Figure. We can write va in 
terms of an anisotropy parameter Ao: 

VA = tan-\^^), (3.1) 

where Ao can be evaluated in terms of Fermi surface parameters. For an ellipsoidal FS one has Ao = y 3 -, but our 
results for the nonlinear currents are not restricted to this case, since they depend only on the properties at the nodes, 
not on the details of the shape of the FS. Then, vj„ = Vt n x n and Un = Vf2 = u/3 = w/4 with the directions given by: 

{ca-Sa), (cA + S A )t - - (ca-Sa), . (cA+S A )t , 

xi = -X3 = -= a -= b, x 2 = —X4, = -= a H -= 6, (3.2) 

V2 v2 v2 v2 

where a, b are unit vectors along the a, b crystallographic axes, ca = cos 1/4 an< i S A = siau A - The calculation is then 
quite analogous to that done in the previous section. We can obtain, from Appendix |X[ j gp j as given in Eq. ( 2.1 5| ) 



with Vf being replaced by Vf n . The linear part is obtained most conveniently in terms of the principal axes of the 
superfluid density tensor, which in this case are the crystallographic directions. If we now define the dimensionlcss 
length Z in terms of A„ (A~ 2 = ^p-iV /««/„), the penetration depth at the nodes, (instead of X a b), v c = and 
Hq = e £ A * f , then by following the steps leading to the derivation of ( 2.2 2j ) and the discussion from Appendix |A], we 



obtain: 



A 2 

dzzUi - -rf Ui + [eu(ip)ul + e 2 i(ip)u a u b + e 3i (ip)ul] = 0, i = a,b (3.3) 



which is formally identical to (2.22) and is solved in the same way. We can express the linear part of the velocity 



components along the principal axes as: 

Li(Z) = at sinh(riZ), i = a,b (3.4a) 

[ir a cosh(r a /4 s ) fin cosh(nZ s ) 

where we have defined = A n /Aj, i = a,b. The nonlinear part satisfies: 

dzzN t - r 2 Y 4 + [eual sinh 2 (r a Z) + e 2l a a a b sinh(r Q Z) smh(r b Z) + e 3i al sinh 2 (r b Z)], (3.5) 
with the solution N^Z) = N hl {Z) + N pi (Z), where: 

N hi (Z) = A u smh(nZ) + A 2i cosh(r l Z), i = a,b (3.6) 

N pi (Z) = cosh(2r a Z) + R, cosh((r a + r b )Z) + S, cosh((r a - r b )Z) + T t cosh(2r b Z) + U t . (3.7) 
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After determining the unknown coefficients from the boundary conditions, as explained in Appendix 0, we obtain: 



N a (Z s ) 



12 



ei 



cos 2 (^-f) 3e 2a cos(^-f)sin(</>-f) 3e 3a sin 2 (^ - f ) , 



r\r h (2r a + r b ) 



r a r 2 (r a + 2r b ) J ' 



(3.8a) 



Nb(Zs) = tttI 



h 2 t 3e lb cos 2 {ip - f ) 3e 2h cosO - f )sin(V> - f ) e 3() sin 2 (-0 - f ) , 



12 r 2 r b (2r a + r fc ) 



r a r 2 (r a + 2r&) 



(3.8b) 



In addition to the nonlinear transverse magneti c mo ment, quad ratic in H a , there is a transverse component linear in 
H a due to the anisotropy of the FS. From Eq. (2.26) and (2.34) we get for the linear part, rrv™ {ip): 



mT(ip) = ^-H a (X a - \ b ) cos 2^. 

47T 



(3.9) 



Experimentally, the different angular and, particularly, field dependence should provide a way to separate the effects 
of m] n from the quantity of interest, the nonlinear m± which reflects the symmetry of the superconducting gap, 



rather than that of the normal state. From Eq. (2.34) and (3.8) we can express the nonlinear part of m±, as in the 
d + s case, in terms of the angle ip: 



SX n H 2 . 1 (c A + s A ) 2 



{c A ~ S A ) 2 



(18 



(c A + S A ) 2 

r a r 2 (r a + 2r b ) 



(CA ~ SaY 



) cos ip sin ip] , ip G [ip±, —ipt + 



(3.10a) 



m±(tp) 



SX n H 2 

'~T7> U~( CA + SA)[COS1p 

iZir Hq 

(ca - sa) 2 



(18 



rln{2r a + r b ) 



■ n l ro (ca-sa) 2 

sin?p -Li- — — - 

V/ 4 L rlr b {2r a +r b ) 



(c A + sa) 



2 r ) cos ip sin ip], 



Ip G [—Ipl + 



(3.10b) 



where tan(^i — j ) = v 6, tan^ + z^), as explained in Appendix B. The effect of FS anisotropy on mj_(ip) is similar 
to that due to the presence of an s-wave gap, discuss ed in the previous section. The period is doubled from the pure 
d-wave and isotropic FS value. The results from Eq. (3.10) can be simply extended to include an admixture of s-wave 
gap and we will defer their discussion for the following subsection. 



2. BSCCO type orthorhombicity 



Here we consider the nonlinear current response for a d-wave OP when there is a FS anisotropy consistent with 
the orthorhombicity of materials in the BSCCO family. The effects of the orthorhombic distortion relevant to the OP 
symmetry are discussed in Refs. El^ and in the Introduction. The nodal lines of the d x 2_ y 2 gap functions are along 
the X, Y axes. These axes are now the principal orthorhombic directions. As a consequence, the Fermi speeds at 
the nodes are nonidentical (vfx V/y) anc * the nonlinear current which is generically oc Vj n will differ along the X, 
Y directions. Hence, some small changes are needed in the procedure used for the YBCO case. Wc defin e th e unit 
of length in terms of A = yj X„x X n y , where X n i, i — X, Y is defined analogous to A„, given above Eq. (3.3). We 
introduce A and Vf = ^vJxvJy, in the previous definitions of the dimcnsionless units, for example, in Z = z/X and 
in H from Eq. J2.21]). We define q, = (A/A,) 1 / 2 , i = X, Y and obtain: 



dzzUi 



■ p 2 [e li (ip)u 2 x + e 2l {ip)u x u Y + e 3l (ip)uy] = 0, 



qju t 



2 2 Af_^£ 



X,Y 



(3.11) 



where the only nonvanishing coefficients e\^,3,i are (cf. Eq. ( [2.17 )), e\x = e 3Y = 1 for ip G f0,7r/2] and — e\x = 
e^Y = 1 for ip G [tt/2, it] . The part of the velocity linear in h, Lj, has the form given by Eq. ( p!4| ) with replaced 
by qi and ip — j by ip. The nonlinear part in h, Ni, from Eq. (3.11) satisfies : 



dzzNi(Z) - <gNi(Z) + a 2 sinh 2 ( % Z) = 0, 



i = X,Y 



(3.12a) 







a x = Px e ixa 2 x , ay = Py^yOy, 



(3.12b) 



It is straightforward, following the YBCO case, to write the complete solution for Ni and obtain for a slab: 

i 2 2 l2 2 

Nx(Z s ) = — ^e lx cos 2 i>, N Y (Z S ) = — ^e 3 ysm 2 V>. (3.13) 
Iz q x iZ q Y 

The corresponding nonlinear transverse magnetic moment is 



SX H 2 v 2 v 2 

m±(ip) — — — — - cos ipsmipleix^- cos?/; — e 3 y^r smV>], (3-14) 
12?r H q x q Y 



where e\x and e 3 y are as given below Eq. ( p^ ). We see that the two terms, which correspond to two "jets" 



of quasiparticle backflow excitations in the nodal directions are of different strengths since px /qx Py / Qy ■ We 
illustrate this effect on the angular dependence of m± in Fig. EL We use here the same conventions as in Fig. |^. 
The solid line represents results obtained for m±(tp) in the isotropic FS case, with its maximumjiorrnalizcd to unity. 
They should be contrasted with the broken line which corresponds to results with Ay /Ax = 1.1 ,E3 normalized by the 
same factor. We have used Ref. ^8| to estimate Xi/X n i ~ V%, and set Vfx/vfy — Ay/Ax = 1.1- Half a period for 
m± is shown, the plot can be extended by odd parity over the whole range, ip € [— 7r/2, 7r/2]. The behavior of m± is 
somewhat similar to that for a d + s gap and isotropic FS (Fig. ||) , but instead of rax vanishing for a field applied 
along the antinodes it vanishes here when the field is applied along the nodal directions. 

B. Anisotropic d+s 



The derivation of (|3.10| ) can be combined with that of ( p. 37 ) to include an s-wave component in the energy gap, i.e. 



A(j+s(0) from Eq. ( |2.10| ), in the presence of FS anisotropy of the YBCO type. Examination of the derivations of the 
equations just mentioned shows that the crucial point about the additional s term is that the shifts in the direction 
of the vj„ from the ±X, ±Y axes, arising from the anisotropy and from the s component must be combined. Thus 
the total shift angle va, introduced in Fig. [|, must include the effect of v as well as that of A 2 ,. For an ellipsoidal FS 



this is achieved by replacing A 2 , in Eq. (3.1) with 



A 2 = Aptan(— + u), (3.15) 

and we note that if v is small and Ag sa 1 then v& rs v\ + v, where v\ is the contribution ar ising solely from FS 
anisotropy. With these generalizations, the results for m±_{ip) in terms of va are still given byEq. ( fjgj ). Some care has 



to be taken with redefining dimensionless quantities. The characteristic velocity is taken to be v c = (A 2 — A 2 ) 1 / 2 jv.t 



In Fig. ^ we show the results from Eq. ( |3.10 ) for the nonlinear transverse magnetic moment m±(ip) (the linear 



part is given by Eq. (3.S)). In the Figure, we give results for Ao = \ a /\ — 1-5 (which is in the correct rang for 
YBCO) and take A 2 = A a A;,. Again, the solid line is our standard reference result: pure <i-wave and isotropic FS, 
normalized to unity. The broken lines are labeled by their va values and are normalized by the same factor. For 
the value of A employed, the range of va that we have chosen to display would correspond, for an ellipsoidal FS, 
to slightly negative values of v. This seems physically appropriate to describe YBCO, where the chains are in the b 
direction and A a > A;,. The departure of mj_(V>) from 7r/2 periodicity can be very significant when the effects of v 
and v A are in the same direction. The precise shape of m±(ip) is the result of the combined effects of FS anisotropy, 
reflected in the penetration depth, and of the s wave contribution. Experimental results for mi(^) would yield va, 
and information about the FS (i.e., A ) would be needed to infer v. 

C. Extended s-wave states 

The results we have discussed so far were obtained for nodes along four directions on the FS. It is straightforward to 
incl ude symmetries of the OP with a larger number of np.des, as we have already mentioned in the derivation leading 



to (2.7). Among the allowed spin singlet pairing statesa of a single CuOi plane with square lattice structure, that 
have not been considered yet, are the A2 g or "g-wave", and the "extended s-wave" states, with eight nodes.Ej The 
gap function for the former case can be written as 

A es (0) = A s + A 9 sin 40, (3.16) 
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where |A S | <C |A g |. The g-wave is the limit A s = 0. We have used an angle <p measured from the g-wave noda l 
direction, along the +X axis (more conventionally, for an angle measured from-,the a-axis, the sine in Eq. 
would have been replaced by a cosine, and reflect the same symmetry properties!!] of A es and the simple s-wave gap 
for a square lattice). The differences between pure g-wetve and extended s-wave, written in this form, are similar to 
those between d and d + s gap functions as explored in Section For example, the nodes are shifted from the g-wave 
gap by an angle ±z/, where we have: 



1 . _ lr A 

i sm ( a^ 



(3.17) 



and v is the angle between the direction of closest node to the +X axis and the +X axis. The gap function at the 
nodes can be approximated by 



4(A 2 - A 2 )V^„ 



(3.18) 



with </>„ taken with the respect to the node n of A es (</>). If we take an isotropic FS, the Fermi velocity at the nodal 
directions is given by Vf n x n with 



Xi= -x 5 = cX - sY, x 2 = -x 6 = c a X + s a Y, 
x 3 — — x-j = sX + cY, £4 = — .t 8 = — s a X + c a Y ', 

where c = cosi^, s = sin^, c a = cos a, s a = sin a, and a = tt/4 + v. 



(3.19) 



After obtaining the nonlinear current by performing the sum in (2.7) over all eight nodes, the solution for the 
velocity field is readily found as in the d + s case, noting that now /1 = 4 and the coefficients 61,2,3 i & r e different. The 
period of m±(ip) becomes it/2. The nonlinear current is of the form ( 2.1 5| ) with 



eix = c -s +c a + s a , e 2X = 2(-cs - c s + c a s a - c a s a ) = 2e 1Y 
e 3 x = cs 2 - c 2 s + c 2 a s a + c a s 2 a = \e 2Y , e 3Y 



. c 3_ s 3 „3 , „3 



(3.20a) 



eix = c J + s 3 + c 3 + 4, e 2X = 2(cs 2 - c 2 s + c 2 a s a - c a s 2 a ) = 2e 1Y , 
Z3X = cs 2 + c 2 s + c 2 a s a + c a s 2 a = 2 e 2V> e 3Y = c 3 - s 3 - c 3 + s 3 , 



TT TT 



With the aid of Eq. (2.36) we can obtain the corresponding m±(ip): 

SX H 2 s c 

mAip) = ° tt~[( c + s )( cs 7=( c - s )) cos3 ip + (s - c)(cs — (c + s)) sin 3 ip 

y ' 48tt Ho y/2 V2 

3 

+ co$ipsmip[((c — s)(l + 3cs) + v2c(l — -(c 2 — s 2 ))) cos ip 

+ ((c + s)(l - 3cs) - V2s(l + ^(c 2 - s 2 ))) sinV]], i/> G [v, | - v] 



(3.20b) 



(3.21a) 



m±(ip) 



SXab H 2 



,, [( c _ s )( cs + —J=( c + s )) cos3 ip + ( c + s)(cs H =(c - a)) sin 3 ip 

TT -Hq v2 v2 



+ cos^sin^[((c + s)(l - 3cs) + V2c(l - -(c 2 - s 2 ))) cos-0 



((c -*)(! + 3cs) + %/2s(l + -(c 2 - s 2 ))) sin V]], 



r 7T TT 

^ G [- - ^ - + H 



(3.21b) 



where in the first definition of i?o, in Eq. (2.21), we have replaced A e ff by (A 2 — A 2 ) 1 / 2 . 



D. Mixed d+g Gap Functions 



Mixed pairing states compatible with the presence of orthorhombic distortion also includeoa d + g-wave states, 
i.e., "s~ + d xy " for the YBCO family and "s~ + d x 2_ y 2 n for BSCCO. We consider here only the first case with an 
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isotropic FS. The other cases can be treated similarly. We again assume that the d- wave character is dominant and 
we denote the gap function we use, with ^> A g : 

A d +g(» = A d sin2^ + A g sin4<^, (3.22) 

where <p is measured from the a axis. At the four nodes we can approximate the absolute value of the gap function, 
needed to calculate j qp , with : 

\&d+ g (4>n)\ » |2(A d ± 2A s )0 n | . (3.23) 

We have again contributions from two inequivalent quasiparticle "jets" , in this case due to the different values of the 
gap near the nodes. The solution for the linear part of the velocity is the same as for the pure ci-wave, discussed in 
the previous Section. If we define itj = Vi/2v c , v c = Ad/vf the nonlinear part of velocity satisfies: 

d ZZ N t {Z) -Ni{Z) +a 2 i smh 2 {Z) = 0, i = X,Y (3.24a) 



'A" 



1 + 2r 



eixa x , 



1 - 2r 



(3.24b) 



where r = A s / Ad, and ax, ay are given by Eq. (2.27). The solution of Eq. (3.24) for a slab is readily obtained as: 



Nx(Zs) 



e-xx 



12 (1 - 2r) 

The corresponding transverse magnetic moment is 



cos 2 V, N Y (Z S ) 



h 2 



e 3Y . 2 , 

sin ip. 



12 (1 + 2r) 



S\Hl . e lx . 

— — cos^sm^ cosV 1 - t 

127T zip 1 + 2r 1 



e3F 
- 2r 



sin tp] 



(3.25) 



(3.26) 



where e±x, egy a re given from Eq. (2.17). This result is formally similar to that for an anisotropic ci-wave for BSCCO, 
recall Eq. (PH) and Fig. |[ with the substitutions Px/Ox = V(l + 2r) and p y /Qy = 1/(1 ~ 2r). For brevity, we 
have not included p lots in this and the previous subsections, but the interested reader can easily generate them from 
the analytic forms ( 3.21 ) and ((3.26 ) . 



IV. NODELESS GAP FUNCTIONS 



In this section we discuss the nonlinear supercurrent response when there are strictly speaking no nodes but there 
are regions ( "quasinodes" ) with a very small. although finite superconducting gap. For an example of a nodeless 
gap function of this type we choose s + idoEil pairing. This is a frequently considered candidate for a gap function 
without nodes, and it clearly illustrates the difference in the nonlinear response when the nodes are absent and only 
small minima in the gap function exist. The range of |A S /A<j| that we consider is restricted to the very small values 
compatible with the bounds set by experiment. For these values, the quasiparticle contribution is dominated by the 
behavior at the "quasinodes" and the general procedures presented in the previous sections can be used. We will 
again concentrate here on the nonlinear magnetic moment, as a probe to identify the symmetry of the gap function. 

IS 

For brevity, we will assume an isotropic FS (equivalent nodes) . Anisotropy can be added using the methodology of 
the previous section. We consider an s + id gap function of the form A s+ id(4>) = A s + iAd sin(2</>), with A s <C A^, 
where is the angle with respect to the +X axis, which is still minimum of |A s +jd(0)|. The phase space available 
for quasiparticle excitations is reduced compared to the ci-wave case so that the critical angle (recall Eq. (]2.8|)) is: 



The behavior of m± for the s + id state is similar to that which would occur for some other gap functions, such as 
"anisotropic s-wave" , proposed in the framework of interlayer tunneling by Chakravarty et al. J2j or for d x 2_ y 2 +id xy , 



_ (v f -vy-Aj 

4A 2 



(4.1) 



There are no nonlinear effects present if v < vt = A s /vf or correspondingly if H a is smaller than a threshold field 



Ht which we shall see is Ht = A s /A<j. We can again perform an approximate integration of Eq. (2.4) and obtain an 
equation for the velocity field (-0 £ [0, ?]): 
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dzz^ -u % + (it? - u T )0{ Ul - u T ) = 0, t = X, Y (4.2) 

where itj = 3-, v c is as defined for the c? wave state, and ut = ott - - The step function arises from the nonlinear 
effects being present only for sufficiently large fields. 

We proceed with the perturbation calculation.-as previously developed, noting that we have to treat carefully the 
nonanalyticity associated with the step functionEJ in Eq. fl4.2p . To do so, we divide the Z > portion of the slab 
(the solution for negative Z follows from the symmetry) into two regions, the first where Ui > ut or correspondingly 
Z G [Zti, Zs] and the second for u, < Ut, that is, for Z S [0, Zti), where the position Zt% is to be determined. The 
solution for the linear part Li(Z) is the same as in the d-wave case. For the nonlinear part Ni(Z) we have from ( [f.2j ): 

d ZZ N i -N i +L 2 i -u T = 0, Ze[Z Ti ,Z s ] i = X,Y (4.3a) 

dzzNi -Ni = 0, Ze [0, Z Ti ) i = X,Y. (4.3b) 
Ni(Z < Zt) can then be expressed in the form of Eq. ( 2.25| ) and for Z < Zt we can write: 



Ni{Z) = A 3l sinh(Z), i = X,Y (4.4) 

where we have used iVj(O) = 0. After finding the coefficients Pi and Ri in the particular solution for Z > Zt, we 
are left with four unknown parameters An, A 2 i, and the matching point Zti- They are determined from Eq. 
( [2.29 ) and by requiring continuity of velocity, magnetic field, and current at Z = Zti- One has the same number of 



conditions and of unknowns. Straightforward algebra yields: 



A 1{ = - -^-[tanh(Z s )cosh(Z Tl )(l - 2 sinh 2 (Z Ti )) - 4sinh(Z s )] (4.5a) 
a 2 

- ( y + 4) tanh(Z s ) cosh(Z Ti ), 



A 2l = -± cosh(Z Tl )(l - 2sinh 2 (Z Tl )) + (y + 4) cosh(Z Tl ), (4.5b) 



A 3l = - -^[tanh(Z s )cosh(Z Ti )(l - 2 sinh 2 (Z Ti )) - 4sinh(Z s ) (4.5c) 

a 2 

+ sinh(Z Ti )(l + 2cosh 2 (Z Ti ))] + (y + u T )[smh(Z Ti ) - tanh(Z s ) cosh(Z Tl )], 

Pi = ~4> Ri = -4~ v r> sinh(Z Ti W) = -^r, i = X,Y (4.5d) 
6 2 Oi(^) 

where in the last expression we have emphasized the angular dependence of the matching point Zxi- From these 
equations we obtain that Ni(Z s = Z^i) — 0. For the slab, when Z s ^> 1, the expression for Ni(Z s ) is simply: 

1 4?/ 3 A 

N X {Z S ) = l-h* cos 2 V> - 4 + ^-^jMh fx-), (4.6a) 

12 3hcosip cosi/'Ad 

N y (Z s ) = [-U 2 sin 2 * - n| + ^_]©(^ - -^V). (4.6b) 
12 sm^A^ 

We see from these expressions that for very small fields h < Jit, with Kt = A s /Ad the nonlinear effects vanish. For 
each jet of quasiparticles to be excited the applied field has to be h > hTi{tp), where hTi(ip) are angle dependent 
threshold fields, the minimum value of which is the overall threshold field hx- If we now compute mj_ as in the 



derivation of Eq. ( 2.36|) , we find 
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S\ h H 2 

fn±(tjj) — — - [—p- cos tp sin tp [cos tpQ(h — hrx) — smijjQ(h — Hty)} 
Yztt Hq 

+ 3(^) 2 iJ [cos^e(/i - h TY ) -ainipQ(h- h TX )} (4.7) 
A H 2 

- 2(-r^) 3 -2-[coti/;e(/i- h TY ) -tsmipQ(h- h T x)]], 
Ad rl a 

wher e Htx and Hty are the threshold fields along the X and Y directions respectively, which can be read off Eq. 



(4.6). It is apparent from this result that in this case m± has a more complicated angular and field dependence than 



in the d-wave gap (the limit ut = ► corresponds to the pure d-wave), or indeed than in any of the other cases 

considered in the previous sections. In contrast to what occurs with gap functions with nodes, where the applied field 
dependence is simply given by an overall factor of H 2 , the angular and applied magnetic field dependences of m± arc 
now coupled, with the coupling depending on the ratio A s /Ad. This coupling ca n b est be described by noticing the 



scaling property of the transverse moment at any angle, which follows from Eq. (4.7) 



m^(Kh,K^-) = K 2 m ± (h,^-) 1 (4.8) 
Ad I\d 

where k is a scaling factor. 

In Fig. |t] we show our results for the angular dependence of m± at two values of the dimensionless field, h — 0.05 
and h = 0.1. We plot m±(tp) in the range ip € [0, 7r/4] and this can be extended from the odd parity, m^(7r/4 + ip) = 
—m±(ir/4 — ip), to the range of its period, tt/2. The solid curve represents the <i-wave result, such that its maximum 
at field h — 0.05 is normalized to unity. As in our discussion of Fig. ^, we refer the reader to previous worko, for the 
unnormalized values of this quantity. The broken lines represent mj_(^) for various A s /Ad values, and attain their 
maxima at different angles tpM- The range of A s / Ad included is greater than that allowed by experiment, so that the 
scaling (with k = 2) given by Eq. ( |4.8] ) can be illustrated. The threshold effect is already evident in this Figure. The 
nonlinear effects are suppressed for very small applied fields, but recover quickly and indeed the maximum value of the 
plotted quantity is slightly enhanced at larger fields. This slight enhancement arises in part from the same reasons as 
in the d + s case, discussed above. Although the "is" component leads to a reduction in the phase space available, this 
reduction is made up for, at larger fields, by the suppression of whichever one of the two "jets" is below the angular 
depe ndent threshold, so that the region where the contributions of the two nodes tend to cancel is diminished. I n E q. 



(4.7), for d wave, we see that for < ip < ir/2 both jets contribute. In that limit (from the first line in Eq. (|4.7| )) 
the contributions of two jets have opposite signs, partially cancelling each other. With the addition of a small s-wave 
component there are two trends: first there is a small reduction in the phase space available for the quasiparticle jets, 
resulting in smaller contributions of individual jets, but also there will be a region where cancelation of the two jets 
will be reduced. 

We show next in Fig. || , the variation with applied field of the maximum value of the transverse moment a as 
a function of angle, for various values of A s /Ad. There are several interesting features in this Figure. The field 
dependence is no longer cx h 2 , except for the A s — > limit (pure d-wave, depicted by the solid line). The broken 
lines (corresponding to several values of A s / Aj_ as in the previous Figure) display the existence of the threshold field 
Iit- The value of hx can be clearly seen from the Figure to be A s /Ad, as given above. There is also, as previously 
discussed, a region of parameter space h, A s /Ad, with a slightly enhanced m±^ max compared to the <i-wave case. 



V. CONCLUSIONS 



We have shown in this paper how, at low temperature, the nonlinear magnetic moment can serve as a high quality 
bulk probe to do "node spectroscopy" in HTSC's, that is, to investigate the position of the nodes (or quasinodes) in 
pairing states with possible mixed symmetries. We have obtained analytic results for mj_(V>) and a variety of pairing 
states. Analytic results are possible because the low temperature nonlinear response due to quasiparticle excitations 
depends on the local properties of the FiLaiid of the energy gap near the nodes. We have examined these effects on 
m± comparing them to previous resultsBta for a pure d-wave and isotropic FS. For gaps with nodes, m± is cx H 2 
and we have concentrated on the angular dependence of m±(jp). 

In sections II and III, we have shown that even a small (real) admixture of an s component to the dominant 
d-w&ve OP, would be detectable in mj_ (?/>). We have found that m±(ip) is closely linked to the direction of v/ at the 
nodes (quasinodes) in the energy gap. Even for an isotropic FS, the directions of the V/ at the nodes are no longer 
mutually orthogonal: the nodes become inequivalent and the period of mi(i/)) doubles to n. Experimentally, as can 
be seen from Fig. [| it may still be convenient to examine the tt/2 and tt Fourier components of m±(ip) to determine 
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the admixture of s-wave. We also have studied the effect of the-a — b plane penetration depth anisotropy on the 
pure d-wave state. If experimental values as reported for YBCCO are used, such an anisotropy significantly alters 
the corresponding d-wave result with an isotropic FS. We have included the influence of orthorhombic distortion on 
mx(ip)t for both the YBCO and BSCCO families of cuprates. These effects can be understood in terms of an angle 
va reflecting the combined effects of an anisotropic FS and the admixture of s component. 

In Section |[y| we have investigated gap functions with quasinodes. The concrete example of an s + id gap clearly 
shows that m± is a very sensitive probe that allows one to distinguish nodes from quasinodes. To do so, one must 
consider both its angular and field dependences. The angular dependence is no longer independent of the applied field, 
and the field dependence is not cx H%. There is a characteristic threshold field Ht = HoA s /Ad required to create 
quasiparticles. No nonlinear effects are present for H a < Ht- These distinct features of m±(ip) strongly suggest that 
they are generic to the other forms of gap function with quasinodes. 

We have focused here on clean superconductors. This is because, as shown in the extensive discussion of Refs. 
[T^,|l8|, impurities are unimportant, for clean available samples, except at very small applied fields. The methods that 
we have described in this paper can be readily extended to other forms of gap functions, such as those that attribute 
an explicit role to the presence of chainst3 in YBCO-like cuprates, or to other combinations not included here. 

Possible future extensions of this work also include incorporating an explicit time dependence of the applied field, 
since the use of more sensitive experimental AC techniques to measure the nonlinear effects is contemplated!^] 
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APPENDIX A: CURRENTS 



We give here the general analytic expression for j 9P (v) at low temperature. For a superconductin g ga p that can be 
approximated at the node n by |A(0„)| ss \fiA e ff^ n (j) n \, n = 1, 2, elementary integration of Eq. (2.7) yields: 



J 9 p( v ) ~ ~2ey2 N fnVfnXnV&eff,n / [<t>\ 



2 V f ^eff,n 



en 



(Al) 



n = 1,2,... 



where for various directions of the applied field, summation will be over different nodes and (f> cn is given by Eq. (2.14). 
For YBCO, including the effect of FS anisotropy, we set A e ff yTl = A e ff (see section III), and obtain that nodes 3, 4 
contribute to j qpa ,b if ip G lipit ^]- A t %j) = ipi only node 3 contributes to the quasiparticle excitations, therefore the 
expressions for jq pa ,b from Eq. (Al) with summation over nodes 2, 3 and over nodes 3, 4 should coincide. Writing 
explicitly jq pa , quadratic in h, in dimensionless form we have: 



e la (i>)L 2 a (Z) + e 3a (^)L 2 b (Z) = e 2b (^)L a (Z)L b (Z), i> = ^ 



which gives (we recall Eq. (3.4)) 



tan(V>i(£) - -) = 



A& ,7T 

-x: tail( 4 



va) 



tanh(r a Z) 
tanh(rbZ) 



(A2) 



(A3) 



Unlike in the case of an isotropic FS, the angular range of the applied field contributing to certain nodes is a function 
of Z. For a slab we can simplify various expressions by taking ipi{Z) ~ ^(-Z^) = and this was done in Eq. (3.10). 
In an analogous way we obtain that = —ipi(Z) as well as the range for ip where other nodes contribute. This 

approximation produces a very small discontinuity in the results at the end points of the intervals in Eq. 2.37. This 
has been interpolated over in Fig. |(| 



15 



APPENDIX B: NONLINEAR VELOCITY FOR AN ANISOTROPIC FS 



After evaluating jg P , as discussed in Appendix^, we ca n ca st Eq. (2.5) in dimensionless form. From the perturb ation 
meth od the nonlinear part in velocity, N, satisfies Eq. (3.5), which is valid for any ijj. The coefficients in Eqs. ( |3.6| ) 
and (3/7) are: 



, a , w ry \ 2r Q sinh(2r Q Z s ) r a + r b sinh((r a + r b )Z s ) 
Mi = -A2it&nh(riZ s ) Pi — — — R4 



Ti cosh^Z,) 
_ r a - r b ^ sinh((r a - r b )Z s ) _ 2r b ^ sinh(2rj,Z s ) 
Ti 1 cosh(riZ s ) ri % cosh(riZ s ) 

A 2i = -Pi - S t -Ti-Ui, i = a,b 



cosh(rjZ s 



(Bla) 



h 2 ei a cos 2 w 
~V 2 6ri cosh 2 (r a Z s )' 



Pi, 



h 2 e lb 



cos 2 U) 



M 2 2r 2 (4r 2 -r 2 )cosh 2 (r Q Z s )' 



(Bib) 



Ra 

R b 



h 2 e 2a 



cos oj sin uj 



H 2 2r a 2r 2 (2r a + n) cosh(r a Z s ) cosh(r h Z s ) 



h 2 e 2b 



cos lu sin lu 



y, 2 2r 2 r b {r a + 2r b ) cosh(r a Z s ) cosh(r b Z s ) ' 



(Blc) 



S a 
S b 



h 2 e 2a 



cos uj sin uj 



fi 2 2r a r 2 (2r a - r b ) cosh(r a Z s ) cosh(r b Z s ) ' 
h 2 e 2b cos w sin uj 

H 2 2rlr b (r a - 2r b ) cosh(r a Z s ) cosh(r b Z s ) ' 



(Bid) 



h 2 e 3a 



sin 2 uj 



fi 2 2r 2 (4r 2 b -r 2 a ) cosh 2 (r b Z s y 



T h = - 



h 2 e 3b sin 2 uj 
M 2 6r 4 cosh 2 (r b Z s 



(Ble) 



h 2 . ei a cos 2 uj 
< \ 



e 3a sm uj 



[i 2 2r 2 a V 2 cosh 2 (r Q Z s ) r 2 cosh 2 \r b Z s ) 



), 



h 2 e\ b cos 2 uj 



e 3a sm a; 



M 2 2r 2 V 2 cosh 2 (r Q Z s ) rg cosh 2 (r b Z s ) " 



(Blf) 



where uj = ?/) — ? and the coefficients e,- 1.2,3 can be obtained as in the isotropic case. 
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FIG. 1. Coordinates and definitions for the d + s OP calculations. The FS and the energy gap are shown schematically. The 
axes a and b are along the crystallographic directions. The orthogonal X and Y axes are alon g the nodal directions of the pure 
d-wave gap. The d + s nodal directions, labeled 1,2,3,4 are shifted by an angle ±v (see Eq. ( |2.1l| )) from their A s = values. 
The applied magnetic field, H a , forms an angle ijj with the +Y direction. 



FIG. 2. Effect of an s component on the angular dependence of the nonlinear transverse magnetic moment. The quantity 
plotted is m±(ip), calculated for a d + s gap and an isotropic FS, at constant field, normalized so that its maximum is unity 
f or th e pure d-wave (solid line). The broken lines are labeled by the corresponding values of the quantity v as defined in Eq. 
( 2.11). A full period is shown. 



FIG. 3. Fourier components for the case shown in Fig 
the 7r (broken line) component, at a constant field, as a function of angle v ( Eq 
7r/2 component is unity for a pure d-wave (v = 0). 



|. The results shown are the 7r/2 Fourier component (solid line) and 
(2.11)). Normalization is taken so that the 



FIG. 4. Definitions of various quantities for an anisotropic FS with YBCO type orthorhombicity. The axes a, b, and X, Y, as 
well as the angles v and ip, are defined as in Fig. [j]. For clarity, only one node (node 2 in the scheme of Fig. |l|) is shown here 
and the OP is not depicted. The angle va (at the node shown) is the angle between +Y and V/2. It has the same magnitude 
at the other three nodes. The vector V is the superfluid velocity. 



FIG. 5. Angular dependence of m±(ip), for a pure d-wave OP and a FS with BSCCO type orthorhombicity . The quantities 
plotted are normalized as in Fig. ^. The solid line represents the isotropic case, and the dashed line the anisotropic case with 
\y/\x = 1.1. 
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FIG. 6. Angular dependence of m±(ip), for ad + s OP and YBCO type orthorhombicity. The solid line is the normalized 
reference result (pure d-wave gap and an isotropic FS). The broken lines, each labeled by the corresponding value of the angle 
va (see text and Fig. ^]) are all obtained for \ a /\b = 1.5. 

FIG. 7. Angular dependence of m±(ip), for various admixtures of s-wave component in an s + id energy gap. Results are 
shown for two values of the dimensionless field, h=0.05 (panel (a)) and h=0.1 (panel (b)). The normalization is taken so that 
the pure d-wave (solid line) maximum is unity at h=0.05. The broken lines, are labeled by the corresponding ratio of A s /Ad- 
The scaling property (with k = 2) from Eq. ([0]), can be verified by comparing the results for h=0.05 and those for h=0.1 
with the appropriate ratios A a /Ad increased by a factor of two. 

FIG. 8. Field dependence of the maximum of m±(ip) for various ratios of A 3 /Ad in an s + id state. The solid line is the 
pure d-wave result normalized so that its value is unity for h=0.05. The broken lines corresponds to the same values of A s /Ad 
as in Fig. 7b. For each value of A a /Ad the threshold field hr = A 3 /Ad can be read off from the graph. Except for A s = 0, 
the lines are not parabolic. 



TABLE I. Spin singlet pairing states in a GuO-z plane. 



Symmetry Class 


Informal name 


Representative State 
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s ( S + ) 


const. 


A 2g 

Big 
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B 2 g 


d xy 


xy 
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Fig. 1: I. Zutic and O.T. Vails, Gap Spectroscopy 
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Fig.2: 1. Zutic and O.T. Vails, Gap Spectroscopy 




Fig. 3: 1. Zutic and O.T. Vails, Gap Spectroscopy 




Fig. 4: I. Zutic and O 
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Fig.5: 1. Zutic and O.T. Vails, Gap Spectroscopy 




Fig.6: 1. Zutic and O.T. Vails, Gap Spectroscopy 




Fig.7a: I. Zutic and O.T. Vails, Gap Spectroscopy 
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Fig. 7b: I. Zutic and O.T. Vails, Gap Spectroscopy 
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Fig. 8: 1. Zutic and O.T. Vails, Gap Spectroscopy 



